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Abstract 



The tt charge asymmetry is measured in events containing a charged lepton (elec- 
tron or muon) and at least four jets, one of which is identified as originating from 
b-quark hadronization. The analyzed dataset corresponds to an integrated luminos- 
ity of 5.0 fb^^ collected with the CMS detector at the LHC. An inclusive and three 
differential measurements of the tt charge asymmetry as a function of rapidity trans- 
verse momentum, and invariant mass of the tt system are presented. The measured 
inclusive tt charge asymmetry is Aq = 0.004 ± 0.010 (stat.) ± 0.011 (syst.). This result 
and the three differential measurements are consistent with zero asymmetry as well 
as with the predictions of the standard model. 
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1 Introduction 

The top quark offers an excellent opportunity to search for departures from the standard model 
(SM) as its large mass makes it unique among all quarks. A possible hint for new physics con- 
tributions showing up in the top-quark sector is the discrepancy of the measured tt forward- 
backward asymmetry with SM expectations, reported by the CDF [IJ and DO [2j Collaborations 
at the Tevatron. These discrepancies are of the order of two standard deviations and even more 
in certain phase space regions. They have generated a large number of theoretical explana- 
tions that attribute them to contributions from physics beyond the standard model (BSM). An 
overview of the variety of theoretical explanations can be found, e.g., in Ref. [3J and references 
therein. 

The production of tt pairs at leading order (LO) is symmetric with respect to the exchange of 
the top quark and antiquark. At higher-order calculations, QCD radiative corrections to the 
qq — > tt process induce an asymmetry in the differential distributions of top quarks and anti- 
quarks. The interference between initial-state and final-state radiation (ISR and FSR) processes 
as well as the interference between the Born and box diagrams generate a correlation between 
the direction of the top-quark momentum and that of the incoming quark, while the direc- 
tion of the top-antiquark momentum is related to that of the incoming antiquark ||4|. While 
these processes induce a forward-backward asymmetry at the Tevatron pp collider, the charge- 
symmetric pp collisions at the Large Hadron Collider (LHC) result in a different effect. At the 
LHC, the larger average momentum fraction of the valence quarks leads to an excess of top 
quarks produced in the forward and backward directions, while the top antiquarks are pro- 
duced more centrally. This makes the difference of the absolute values of the rapidities of top 
quark and antiquark, A\y\ = \yt\ — \yt\, & suitable observable to measure this tt charge asym- 
metry. The rapidity is defined asy = ^ where E denotes the particle energy and pz 
its momentum component along the beam direction. For a given sensitive variable, the charge 
asymmetry is defined as: 



N+-N- 

where and represent the number of events with positive and negative values in the 
sensitive variable, respectively. 

Recently, the Compact Muon Solenoid (CMS) and ATLAS Collaborations have published first 
measurements of the charge asymmetry at the LHC and found respectively Aq = —0.013 ± 
0.028 (stat.)+°;°^^ (syst.) [5J and Ac = -0.018 ± 0.028 (stat.) ±0.023 (syst.) IH, consistent 

with the SM prediction of A*"""^ = 0.0115 ± 0.0006 [7J. Although these results do not have 
the precision to establish a non-zero asymmetry at the LHC, they seem to disfavor large pos- 
itive deviations from the SM prediction. The potential disagreement between the Tevatron 
and LHC results might be due to BSM contributions having different effects on the Tevatron 
forward-backward asymmetry and the LHC charge asymmetry [8|. But even at the Tevatron 
the situation is inconclusive as the CDF Collaboration measures a strong dependence of the 
forward-backward asymmetry on the invariant mass of the tt system, while DO results do not 
show such an effect. To shed light on this question, it is crucial to not only measure the in- 
clusive asymmetry but to also measure Ac as a function of suitable variables enhancing the tt 
charge asymmetry in certain regions. 

In this Letter, we report on updated and further developed measurements of Ac, adopting 
the event selection, background estimation, and reconstruction of the tt system from Ref. ||5|. 
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3 Data and simulation 



We present an inclusive measurement and three differential measurements of the tt charge 
asymmetry as a function of the rapidity, transverse momentum, and invariant mass of the tt 
system, using the full 2011 dataset. Each of these three variables is sensitive to a certain aspect 
of the tt charge asymmetry. 

The rapidity of the tt system in the laboratory frame, \y^^\, is sensitive to the ratio of the contri- 
butions from the qq and gg initial states to tt production. The charge-symmetric gluon fusion 
process is dominant in the central region, while tt production through qq annihilation mostly 
produces events with the tt pair at larger rapidities, which implies an enhancement of the 
charge asymmetry with increasing [7]. 

The transverse momentum of the tt pair in the laboratory frame, Pj, is sensitive to the ratio of 
the positive and negative contributions to the overall asymmetry. The interference between the 
Born and the box diagrams leads to a positive contribution, while the interference between ISR 
and FSR results in a negative contribution. The presence of additional hard radiation implies 
on average a higher transverse momentum of the tt system. Consequently, in events with large 
values of pj, the negative contribution from the ISR-FSR interference is enhanced [7J. 

The charge asymmetry is expected to depend on the invariant mass of the tt system, m^^, since 
the contribution of the qq initial state processes is enhanced for larger values of m^j. This ob- 
servable is also sensitive to new physics contributions; potential new heavy particles could be 
exchanged between initial quarks and antiquarks and contribute to the tt production (see e.g., 
Ref. |9| and references therein). The amplitudes associated with these new contributions would 
interfere with those of the SM processes, leading to an effect on the tt charge asymmetry, which 
increases as a function of the invariant mass of the tt system. 

2 The CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid, of 6 m internal diam- 
eter, providing a magnetic field of 3.8 T. Within the field volume are a silicon pixel and strip 
tracker, a crystal electromagnetic calorimeter (ECAL) and a brass/ scintillator hadron calorime- 
ter. CMS uses a right-handed coordinate system, with the origin at the nominal interaction 
point, the x axis pointing to the centre of the LHC ring, the y axis pointing up (perpendicular 
to the LHC plane), and the z axis along the counterclockwise beam direction. The polar an- 
gle 9 is measured from the positive z axis and the azimuthal angle (p is measured in the x-y 
plane. The pseudorapidity is defined as t] = — ln(tan0/2). The inner tracker measures trajec- 
tories of charged particles within the pseudorapidity range < 2.5, while the calorimeters 
provide coverage up to \tj\ = 3.0. The ECAL has an energy resolution of 3% or better for the 
range of electron energies relevant for this analysis. Muons are measured in the pseudorapid- 
ity range |?/| < 2.4, with gas-ionization detectors embedded in the steel return yoke. Matching 
the muons to the tracks measured in the silicon tracker results in a transverse momentum res- 
olution between 1 and 5 %, for pj values up to 1 TeV/ c. Extensive forward calorimetry comple- 
ments the coverage provided by the barrel and endcap detectors. A more detailed description 
of CMS can be found in Ref. LIOJ. 

3 Data and simulation 

The measurements reported in this Letter are based on data taken with the CMS detector at a 
centre-of-mass energy of 7 TeV, corresponding to an integrated luminosity of 5.0 fb^\ To trans- 
late the distributions measured with reconstructed objects to distributions for the underlying 
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quarks, samples of simulated events are used. Top-quark pair events are generated with two 
different generators, either with MadGraph version 5.1.3 |(TT| or with the next-to-leading or- 
der (NLO) generator POWHEG [12J. For both samples the parton shower is simulated using 
PYTHIA version 6.4.24 |13|, and the MLM parton shower/ matrix element matching |IT4| in case 
of MadGraph. Also the t and tW channels of electroweak production of single top quarks 
are simulated using POWHEG. The production of electroweak vector bosons in association with 
jets (W+jets and Z+jets) is simulated using the same combination of MadGraph and PYTHIA 
as for the tt signal. All samples are generated using the PYTHIA Z2 Monte Carlo tune | [T5| to 
model the underlying event. Additional proton-proton interactions (pileup) are overlaid on 
the simulated events as observed in the analyzed data. 

4 Event selection and estimation of background 

The analysis uses tt events where one of the W bosons from the decay of a top-quark pair 
subsequently decays into a muon or electron and the corresponding neutrino, and the other 
W boson decays into a pair of quarks originating jets. We therefore select events containing one 
electron or muon and four or more jets, at least one of which is identified as originating from 
b-quark hadronization. For the reconstruction of electrons, muons, jets, and any imbalance 
in transverse momentum due to the undetected neutrino, £^^*^, we use a particle-flow (PF) 
algorithm |16|. Electron (muon) candidates are required to have pj > 30(20) GeV/c and be 
within \ri\ < 2.5(2.1), while jets are required to have pj > 30 GeV/c and |?/| < 2.4. More details 
on the selection criteria applied to the events can be found in Ref. ||5|. 

In total, 57 697 events are selected, 24 705 events in the electron +jets channel and 32 992 events 
in the muon+jets channel. About 20% of these events are expected to come from background 
processes like W+jets and Z+jets production, the production of single top quarks, and multijet 
production. For the estimation of the background contributions we make use of the discrimi- 
nating power of E^^^^ and of M3, the invariant mass of the three jets with the largest vectorially 
summed transverse momentum. For each of the two lepton channels, these two distributions 
are fitted with a binned maximum-likelihood fit. For the tt signal and the W+jets, Z+jets, and 
single-top-quark background processes, the respective simulated samples are used to model 
the shapes of the E^^^ and M3 distributions, while an approach based on data from sideband 
regions featuring non-isolated leptons is used for the multijet background. Gaussian rate con- 
straints are introduced into the likelihood function for the Z+jets and single-top-quark pro- 
cesses according to the respective NLO cross sections, while the rates of all other processes are 
free parameters of the fit. A more detailed description of the fitting procedure can be found in 
Ref. la. 

The resulting rates for the different processes can be correlated with each other, which has to be 
propagated to the calculation of the statistical uncertainty of the measured tt charge asymmetry. 
The largest correlations are found between the rates of the Z+jets and multijet backgrounds 
(-20%) and between the rates for the W++jets and W"+jets backgrounds (+12%). All other 
correlations among the fit parameters are found to be small. Table [l] summarizes the results of 
the fits, along with their uncertainties. Figure [ijshows the measured £^'®* and M3 distributions, 
with the individual simulated contributions normalized to the results from the fit. 

5 Measurement of the tt charge asymmetry 

The measurement of the tt charge asymmetry is based on the fully reconstructed four-momenta 
of the top quarks and antiquarks in each event. We reconstruct the leptonically decaying W bo- 
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5 Measurement of the tt charge asymmetry 




Ef"" [GeV] M3 [GeV/c'] 



Figure 1: Comparison of the combined lepton+jets data with simulated contributions for the 
distributions in E^^^ (left) and M3 (right). The simulated signal and background contributions 
are normalized to the results of the fits given in Table [Tj 

Table 1: Results for the numbers of events for background (BG) and tt contributions from fits 
to data, along with their uncertainties. The quoted uncertainties are of statistical nature with 
the exception of the uncertainties on the numbers for the single top and Z+jets backgrounds, 
which reflect the widths of the Gaussian rate constraints used in the likelihood fit (see text). 



Process 


Electron+jets 


Muon±jets 


Total 


Single top {t + tW) 


1113 ±338 


1418 ± 505 


2532 ± 608 


W++jets 


1818 ± 227 


1807 ± 290 


3625 ± 369 


W"+jets 


1454 ± 224 


1320 ± 275 


2773 ± 355 


Z+jets 


535 ± 153 


600 ± 170 


1135 ±229 


Multijet 


1142 ±227 


863 ± 209 


2005 ± 308 


Total BG 


6062 ± 540 


6008 ± 698 


12070 ± 882 


tt 


18634 ± 390 


26976 ± 468 


45610 ± 609 


Observed data 


24705 


32992 


57697 



son from the measured charged lepton and E^^^^, and associate the measured jets in the event 
with the quarks in the tt decay chain. The reconstruction procedure is described in detail in 
Ref. [5J. 

The reconstructed top quark and antiquark four-momenta are used to obtain the inclusive (see 
Fig. [2]) and differential distributions of A\y\ and the charge asymmetry is calculated from the 
number of entries with A\y\ > and A\y\ < 0. In case of the differential measurements, the 
asymmetries are calculated separately for the different bins in the kinematic variable y„ where 
Vi is either pj, or m^^. To allow for a comparison of the resulting asymmetry and the 
predictions from theory, the reconstructed distributions of A\y\ and the three kinematic vari- 
ables have to be corrected for background contributions, reconstruction effects, and selection 
efficiencies. 

In the first correction step, the distributions of background processes are normalized to the 
estimated rates (see Table [l]l and subtracted from the data, assuming Gaussian uncertainties 
on the background rates as well as on statistical fluctuations in the background templates. The 
correlations among the individual background rates are taken into account. 

The background-subtracted distributions are transformed from the reconstruction level to the 
particle level after event selection, and from there to the particle level before event selection. 
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Figure 2: Comparison of the combined lepton+jets data with simulated contributions for the 
distributions in A|i/|. The simulated signal and background contributions are normalized to 
the results of the fits given in Table [T| 



The corrections are achieved by applying a regularized unfolding procedure to the data flT) 
through a generalized matrix-inversion method. In this method, the perturbing effects are de- 
scribed by a smearing matrix S that translates the true spectrum x into the measured spectrum 
w = Sx. As reconstruction and selection effects factorize, the smearing matrix S can be con- 
structed as the product of a migration matrix and a diagonal matrix with the efficiencies for 
each of the bins on the diagonal, and all other elements set to zero. The unfolding procedure 
used for the inclusive measurement — described in detail in Ref. 151 — can be generalized to 
deal also with two-dimensional distributions. 

The binning choice for a two-dimensional unfolding procedure has to fulfill some requirements 
in order to stabilize the unfolding procedure and to avoid a loss of resolution. For the applied 
unfolding procedure it is advised to use twice as many bins for the reconstructed spectra as for 
the unfolded spectra [17]. We use 16 (8) bins for the reconstructed (unfolded) A\y\ distribution 
and 6 (3) bins in the reconstructed (unfolded) Vj distributions. Furthermore it is desirable that 
the number of entries in each bin of the reconstructed distributions as well as in the unfolded 
distributions be approximately equal. The ranges for the bins in the unfolded kinematic vari- 
ables are [0 - 0.41; 0.41 - 0.90; 0.90 - oo] for |i/tt|, [ - 23; 23 - 58; 58 - oo] for p^* in GeV/c, and [0 - 
420; 420 - 512; 512 - oo] for JWjj in GeV/c^. The binning choice for A |i/| is different in each bin of 
Vi, resulting in different amounts of vertical overlap between horizontally neighbouring bins 
in the two-dimensional distributions (for illustration see the binning in Fig.|3| lower right). For 
the regularization of these distributions used in the differential measurements all combinations 
of neighbouring bins are considered. Due to the partial vertical overlap of horizontally neigh- 
bouring bins for a given central bin there are up to four possible combinations, each weighted 
with a factor considering the amount of vertical overlap. 

We use separate migration matrices and selection efficiencies for the inclusive measurement 
and the three differential measurements, obtained from tt events simulated with POWHEG. 
Figure |3] shows the migration matrices for the inclusive measurement and for the differential 
measurement in m^l, as an example for the three migration matrices for the differential mea- 
surements. While for the inclusive measurement the migration matrix describes the migration 
of selected events from true values of A|i/| to different reconstructed values, for the migration 
matrices of the differential measurements not only the migration between bins of A |i/ 1 has to be 
taken into account, but also the migration between bins of Vi. The migration matrices for the 
differential measurements feature on large scale a grid of 6 x 3 bins in Vj with each of these bins 
hosting a 16 X 8 migration matrix describing the migration between different A\y \ values. The 
values of A I y I and V, affect the probability for an event to survive the event selection criteria. 
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5 Measurement of the tt charge asymmetry 



The selection efficiencies as a function of A|i/| for the inclusive measurement and as a function 
of A I y I and m^^. for the differential measurement in m^j are depicted in Fig. |3] The nearly sym- 
metric shapes of the efficiency distributions imply that the effect of the event selection criteria 
on the tt charge asymmetry is small. 




Figure 3: Migration matrix (upper row) between the generated and the reconstructed values 
after the event selection for A\y\ (left) and for the measurement differential in m^i (right). Selec- 
tion efficiency (lower row) as a function of generated A|i/|, defined with respect to inclusive tt 
production (left) and for the measurement differential in m^j (right). 

The performance of the unfolding algorithm is tested in sets of pseudoexperiments, each of 
which provides a randomly generated sample distribution. For each pseudoexperiment the 
number of events from each contributing process is determined from a Poisson distribution 
around the mean of a Gaussian distribution centred around the measured event rate given 
in Table [ij with a width corresponding to the respective uncertainty. We randomly draw the 
resulting number of events for each process from the respective simulated sample to generate 
distributions for each pseudoexperiment. Each distribution is then subjected to the background 
subtraction and unfolding procedure described above. 

The average asymmetries from 50 000 pseudoexperiments for the inclusive as well as for the 
differential measurements agree well with the true asymmetries in the sample used to model 
the signal component and the pull distributions agree with expectations, indicating that the 
treatment of uncertainties is consistent with Gaussian behavior. To test the unfolding proce- 
dure for different asymmetries, we reweight the events of the default tt sample according to 
their A\y\ value with a factor if = /c ■ A|i/| + 1, to artificially introduce asymmetries between 
—0.2 and +0.2, and then perform pseudoexperiments for each of the reweighted distributions. 
For the differential measurements this test is performed in each of the three bins of F, separately. 
We find a linear dependence of the ensemble mean on the input value for the asymmetry. In 
addition to this global reweighting of events, one can define the reweighting factor w asa func- 
tion of one of the kinematic variables, w = k{Vi) ■ A\y \ + 1. Four scenarios with k rising or 
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falling linearly with V\ and one scenario in which k rises quadratically are tested, generating 
asymmetries between —0.1 and +0.1. The effect of this reweighting dependent on is tested 
in all three possibilities to measure Aq as a function of Vj. These scenarios serve as tests of the 
model-independence of the unfolding procedure, and observed deviations from the expecta- 
tions are considered for the estimate of the systematic uncertainties of the measurement. 

6 Estimation of systematic uncertainties 

Systematic uncertainties with an impact on the differential selection efficiency, on the recon- 
structed top-quark momenta, or on the background rates can bias the results. To evaluate each 
source of systematic uncertainty, we repeat the background estimation and the measurement 
of Aq using modified simulated samples. The expected systematic uncertainty for each source 
is taken to be the shift in the values of the corrected asymmetry between the default measure- 
ment and the one using the modified templates. The systematic uncertainties can be divided 
into three different categories: experimental sources, uncertainties in the modeling of the signal 
and background processes, and uncertainties due to the applied unfolding procedure. 

The following experimental sources of systematic uncertainties are evaluated: variations in the 
jet energy scale (JES), jet energy resolution (JER), b-tagging efficiency, and the lepton selection 
efficiency. In order to derive the modified templates the corrections on JES and JER for simu- 
lated events are changed by ±1 standard deviations of their ?/- and px^dependent uncertainties. 
The overall scaling factor of the b-tagging efficiency does not affect the measurement, only ?/- 
dependent variations could in principle change the results. We therefore reweight events with 
b-tagged jets in the central region and in the forward regions maximally different within the 
b-tagging efficiency uncertainty. The effect on the Aq measurement is found to be negligible. 
In a similar manner, we vary the scale factor for the lepton selection efficiency within its un- 
certainties, this time with maximally different weighting factors for positively and negatively 
charged leptons, as a possible difference could lead to artificial asymmetries. 

Regarding the simulation of signal and background processes, several sources of systematic 
uncertainties are evaluated. The uncertainty associated with the choice of the event generator 
used for modeling the tt signal is estimated by using simulated events generated with Mad- 
Graph instead of POWHEG for the determination of the smearing matrix. In addition, a signal 
sample with a different hadronization and shower modeling has been used (mc@NLO 4.0 [181 
interfaced to HERWIG ||T9| ) to estimate the systematic uncertainty due to this part of the event 
generation. The effects of variations in the factorization and renormalization scales (Q^) are 
estimated for W+jets and tt events. For this purpose the strong coupling constant as and the 
parton distribution functions (PDF) are recalculated for each event for the varied scale — 
either multiplied with a factor of 4 or 0.25. The scale is varied independently for W+jets and 
tt processes and the estimated uncertainties have been added in quadrature to obtain the result- 
ing systematic uncertainty on the measurement. The systematic uncertainties on the measured 
asymmetry from the choice of PDFs for the colliding protons used in the simulated events are 
estimated using the CTEQ6.6 |20| PDF set and the LHAPDF ||2ll package. In addition, the effect 
of variations in the frequency of occurrence of pileup events, overlaid on the simulated signal 
and background events, is estimated. 

As the first step of the correction of the measured distributions is the subtraction of the back- 
ground, the measurement is sensitive to the asymmetries present in the background model and 
we therefore evaluate the influence of possible mismodeling of the two backgrounds which 
are most sensitive to mismodeling. The rates for positively charged and negatively charged 
W bosons are asymmetric and since the distributions of the two processes are slightly different. 
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6 Estimation of systematic uncertainties 



a mismodeling could artificially produce asymmetries. To estimate the effects from possible 
mismodeling of the W+jets background component, the templates for W+ and processes 
are interchanged. The other background process that can show artificial asymmetries is the 
multijet background. This is the case if the rates for positively and negatively charged leptons 
differ in this sample. The multijet background is modeled using events from a sideband re- 
gion, defined by inverting the requirements on the isolation of the charged lepton candidates. 
In these events the lepton rapidity is on average larger than the jet rapidities. As a result, the 
reconstructed leptonically decaying top quark candidates have on average a larger absolute 
value of the rapidity than the hadronically decaying top quark candidates, which in the end 
leads to different mean values of A|i/| for events with positively and negatively charged lep- 
tons, respectively. To account for this effect, we invert the sign of A|i/| for each event and use 
this altered template to model the multijet background. 

Table 2: Systematic uncertainties for the inclusive measurement of Aq- 



Systematic uncertainty 


Shift (±) in inclusive Aq 


JES 


0.003 


JER 


0.002 


Lepton ID/ sel. efficiency 


0.006 


Generator 


0.001 


Hadronization 


0.001 


scale 


0.002 


PDF 


0.002 


Pileup 


< 0.001 


W+jets 


0.004 


Multijet 


0.001 


Migration matrix 


0.002 


Model dependence 


0.007 


Total 


0.011 



The third category of systematic uncertainties deals with the impact of the limited number of 
simulated events and possible violations of the assumption that the applied unfolding pro- 
cedure is model-independent. The impact of statistical uncertainties of the entries in the mi- 
gration matrices is evaluated by repeating the measurement with altered migration matrices, 
where each element is varied within its statistical uncertainties. In addition to these uncer- 
tainties, we estimate the influence of possible dependencies of the asymmetry on one of the 
three kinematic variables Vj ("model-dependence"). We perform pseudoexperiments with 
reweighted simulated signal samples and evaluate the differences between true and measured 
asymmetries in various reweighting scenarios, as described above. We take the average of the 
observed deviations and assign it as systematic uncertainty. 

The contributions of the different sources of systematic uncertainties to the total uncertainty 
of the inclusive measurement are summarized in Table |2] The total systematic uncertainty is 
smaller than the one obtained in Ref. [5|. The two main changes in the evaluation of systematic 
uncertainties with respect to Ref. [5|, which account for this difference, are discussed below. 
Variations in the threshold for the matching of matrix elements and parton shower evolution 
for the simulation of the tt signal [14], causing the largest contribution to the total systematic 
uncertainty in the previous measurement, have no impact on the present measurement due to 
the usage of the NLO event generator powheg for modeling the tt signal. Furthermore, we do 
not quote a separate uncertainty due to variations in the amount of ISR and FSR on the mea- 
surement, as this contribution is covered by the uncertainties due to the choice of the scale 
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and the model-dependence systematic. The probability for additional radiation increases with 
decreasing and vice versa. Due to the strong correlation between the amount of additional 
radiation and the transverse momentum of the tt system, the variation of the generated asym- 
metry as a function of pj, as done in the estimation of the model dependence uncertainty, is 
also suited to estimate the effects of variations in the amount of ISR/FSR on the measurement. 

The systematic uncertainties on the differential measurements are included in the error bars 
of the corrected differential distributions (see Fig. |4j|. Depending on Vj and the actual bin, the 
contributions from the different sources vary. The largest contributions arise from variations 
in the JES and the conservatively estimated uncertainties due to lepton selection efficiency and 
model-dependence as well as the statistical fluctuations of the migration matrix. The generator 
and hadronization uncertainty play a significant role for the measurements differential in tn^i 
and Pj. The modeling of the major background, the W+jets process, is significant in the third 
bin of JWjj and |. 

7 Results 

The unfolded A\y\ distribution, shown in Fig. [Ij is used to calculate the corrected inclusive 
asymmetry: 

Ac = 0.004 ± 0.010 (stat.) ± 0.011 (syst.) . (2) 

Table |3] gives the values of the measured inclusive asymmetry at the different stages of the 
analysis. 

Table 3: The measured inclusive asymmetry at the different stages of the analysis and the 
corresponding theoretical prediction from the SM. 



Uncorrected 


0.003 ±0.004 (stat.) 


BG-subtracted 


0.002 ± 0.005 (stat.) ± 0.003 (syst.) 


Final corrected 


0.004 ±0.010 (stat.) ±0.011 (syst.) 


Theoretical prediction (SM) 


0.0115 ±0.0006 



The results of the three differential measurements can be found in Table |4] and Fig. |4j The 
measured values are compared to the SM predictions — based on the calculation of Ref. ||7| — 
and as an illustrative example to the predictions from a BSM model that introduces an anoma- 
lous effective axial-vector coupling to the gluon at the one-loop level Il22ll23l| . The gluon-quark 
vertex is treated in the approximation of an effective field theory with an order of 1 TeV scale 
for new physics contributions. This is a model that can explain the strong dependence of the 
forward-backward asymmetry on m^j as seen by CDF. As the theoretical predictions are nor- 
malized to the leading-order cross section and pj is zero at LO, no theoretical predictions are 
available for this differential measurement. Instead, we compare the measured asymmetries 
with the predictions obtained from POWHEG simulation. Within the uncertainties the data do 
not show any significant asymmetry and all measured values are consistent with a null asym- 
metry as well as with the SM predicted values. 
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8 Summary 



Table 4: The corrected asymmetry values in three bins of the kinematic variables |i/tt|/ Pxtf/ ^rid 
nif-i with statistical and systematic uncertainties, along with the SM predictions (in case of pj 
we compare to the values obtained from POWHEG simulation). 



Kinematic variable 


Ac in bin 1 


Ac in bin 2 


Ac in bin 3 


lytti 

|i/tt|(SMpred.) 


0.029 ± 0.021 ± 0.010 
0.0030 ± 0.0002 


-0.016 ±0.015 ±0.010 
0.0086 ± 0.0004 


0.001 ± 0.026 ± 0.022 
0.0235 ± 0.0010 


Pj (simulation) 


0.037 ± 0.025 ± 0.022 
0.0185 ± 0.0004 


0.014 ±0.014 ±0.012 
0.0022 ± 0.0004 


-0.030 ± 0.021 ± 0.019 
0.0006 ± 0.0004 


m^^ (SM pred.) 


-0.051 ± 0.027 ± 0.021 
0.0077 ± 0.0003 


0.017 ±0.017 ±0.014 
0.0112 ± 0.0004 


0.019 ±0.017 ±0.023 
0.0157 ± 0.0006 
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Figure 4: Unfolded inclusive A\y\ distribution (upper left), corrected asymmetry as a function 
of |i/t{| (upper right), pj (lower left), and m^l (lower right). The measured values are compared 
to NLO calculations for the SM — based on the calculations of Ref . [7J — and to the predictions 
of a model featuring an effective axial-vector coupling of the gluon (EAG) |23|. The error bars 
on the differential asymmetry values indicate the statistical and total uncertainties, determined 
by adding statistical and systematic uncertainties in quadrature. The shaded areas indicate the 
theoretical uncertainties on the NLO calculations. 



8 Summary 

An inclusive and three differential measurements of the charge asymmetry in tt production 
at the LHC have been presented. Events with top-quark pairs decaying in the electron±jets 
and muon±jets channels were selected and a full tt event reconstruction was performed to 
determine the four-momenta of the top quarks and antiquarks. The observed distributions 
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were then corrected for acceptance and reconstruction effects. The measured values constitute 
the most precise determination of the tt charge asymmetry at the LHC to date. Although these 
measurements do not allow distinguishing a zero asymmetry from the values predicted in 
the standard model, the reported results nonetheless indicate that LHC data disfavor large 
deviations from the SM predictions. 
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